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Quantum initial state estimation through entanglement and continuous measurement is intro- 
duced. This paper provides a unified formulation of classical and quantum smoothing and shows 
a smoothing uncertainty relation. As an example, a communication between two parties via a two 
mode squeezed state is shown. 
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I. INTRODUCTION 



Quantum filtering theory is an important basis of 
quantum feedback control if control inputs are designed 
from measurement outcomes. A wide range of applica- 
tions of quantum filtering can be found in recent theoret- 
ical and experimental works such as spin squeezing [l|, Q , 
single photon state [|| and superposition state generation 
[3] . The idea of filtering was initiated by Wiener as a part 
of signal processing Q ■ The Wiener filter gives an esti- 
mate if a given signal is detected by noisy measurement. 
While Wiener's work was focused only on a measurement 
process, Kalman considered estimation for a case where 
a system consists of two parts @: One is a dynamical 
part which produces a signal according to the internal 
dynamics of the system, and the other is a measurement 
process of the signal produced from the dynamical part. 
Then, the estimates of system observables are given by 
the Kalman filter. His basic formulation of the system is 
still useful for quantum systems [7|, |8[ . Quantum filtering 
can be thought of as an extension of his formulation and 
given in many different contexts 0, [13, 03 ■ O n the other 
hand, in quantum optics, filtering was developed as an 
expression of a posteriori density matrix under continu- 
ous measurement in terms of quantum trajectory [l2| . 

These quantum filtering techniques are used to obtain 
the expectation of system observables or a density matrix 
at time t conditioned on noisy measurement outcomes up 
to t. However, classical estimation theory includes more 
general situations. In the classical case, given measure- 
ment data up to t, we can calculate a conditional ex- 
pectation of observables at any time. In particular, the 
estimation of past observables at to < t is called smooth- 
ing. Unfortunately, smoothing is not always possible 
in the quantum case. This is because past observables 
are not compatible with present measurement outcomes. 
As a result, the conditional expectation of the past ob- 
servables under continuous measurement is not generally 
well-defined. 

Naturally, however, there will be a situation where we 
want to estimate the initial state of a quantum system 
using continuous measurement. To avoid the compat- 
ibility problem, we have to make a measurement on a 
different system from the one that we want to estimate. 
In this case, there should be a correlation between these 



two systems. In other words, the estimation of the past 
observables is possible only through entanglement. 

For example, this requirement is satisfied by a double 
well Bose-Einstein Condensate (BEC) system [13j ]. Due 
to the symmetry of the double well structure and the 
property of bosons, the BEC state is described by an 
entangled state between the two wells. Then, we want to 
estimate an observable such as the number of particles 
in one well by injecting an optical field and continuously 
detecting the output field. In general, the number of 
particles varies in time [3]. And also, the optical field 
scatters the BEC particles and reduce the number [l5l ]. 
However, due to the initial entanglement between the two 
wells, we can know the number of particles by estimating 
the state of the other well from measurement outcomes 
of the continuous detection. 

Another example is the estimation of past canonical 
observables. Initially two parties, Alice and Bob, share 
a symmetrically prepared entangled state such as a two 
mode squeezed state. Alice cannot estimate her own ini- 
tial state using continuous measurement because of the 
incompatibility between her measurement outcomes and 
the past observables. Instead, she estimates Bob's initial 
state. This procedure can be used as a communication 
between them if Bob's local operation is encoded in the 
shared entanglement. After recording measurement out- 
comes for a certain period of time, she can obtain a pre- 
cise estimate. This can also be thought of as a smoothing 
problem. 

This paper introduces a general formulation of quan- 
tum smoothing. As an example, we consider quantum 
linear systems with Gaussian states and give a detailed 
analysis for a case where the initial state is described 
by a two-mode squeezed state. Then, it will be shown 
that the two past canonical observables can be simulta- 
neously estimated with arbitrary accuracy through quan- 
tum smoothing if the initially shared state is perfect. We 
also give an information theoretical analysis of quantum 
smoothing which clearly shows a difference between a fil- 
ter and smoother. 

We first review classical filtering and smoothing in 
Sec|Hl and then introduce the quantum case in Sec lIIII 
To consider an example of smoothing, Sec lIVI introduces 
a formulation of quantum linear systems and a detailed 
analysis is given in Sec|V] Some singular situation is dis- 
cussed in Scc lVTI 
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II. CLASSICAL CASE REVISITED 

We first introduce classical filtering and smoothing. 
This will help us to understand a quantum analog of 
smoothing because quantum filtering and smoothing can 
be derived along the same line as the classical case. 



A. Conditional expectation 



Let us consider a nonlinear classical system 



dx t = a{x t )dt + B{xt)dwt, 
dm t = c(x t )dt + Ddw t , 



(la) 
(lb) 



where the vector xt represents the state of the system, Wt 
is the input vector and rrit is the measurement outcome 
vector, wt is a classical Wiener process satisfying Ito rule 



dw t dw t 



Idt, 



(2) 



where I is the identity matrix. The first equation de- 
scribes the dynamics of the system and the second one is 
the output process. 

Denoted by P the probability distribution of the 
Wiener process w t - Our purpose here is to calculate the 
conditional expectation of an arbitrary function f(xo, Xt) 



n t [f(x ,Xt)] =Ep[f(x ,x t )\Mt] 



(3) 



where Mt is a measurement space (a closed space gener- 
ated by linear combinations of a constant and {m s |0 < 
s < t}.) This should be normalized as 7r t [l] = 1. 

It follows from Girsanov theorem that if wt is Brownian 
under P, then m t is Brownian under Q defined by 

dP = z t dQ, (4) 

where Zt is a Radon-Nikodym derivative given by 



zt = exp 



1 



^{DD^DdWs + -cHDD^^cds 



(5) 



Thus, by Bayes rule, the conditional expectation is 
rewritten as 



E P [f(x ,x t )\M t 



EQ[z t f(xQ,x t )\Mt 
B Q [z t \M t ] 



(6) 



Ignoring normalization, we can define the conditional ex- 
pectation of / as 



TTt[f{x ,x t )} ~ EQ[ztf(xo,xt)\M t ], 



(7) 



which indicates that the infinitesimal evolution of the 
conditional expectation d7r t [/] is obtained by calculating 
d(z t f). By definition, we have 



dz t = c'(DD') 1 z t dm t , 
df = ?fdt + (Vf)bdw t , 



(8) 
(9) 



where V = d/dx t and 
d_ 

at 



'j 



J2 a * v * + ? J2 / '"- / '// v < v ,'- ( io ) 



hi 



As a result, we obtain 



d(fzt) =(?f)ztdt + [(Vf)BD* + fJ]{Drf)- l z t dm t 

(11) 

The infinitesimal evolution of the conditional expectation 
is then given by 

=n t [f}+TT t [3f]dt 

+ Tr t [(Vf)BD* + fc^DD^dmt, (12) 

where iit+dt denotes an unnormalized conditional expec- 
tation at time t + dt. 

Normalization requires the evolution of / = 1, which 
is given as 



n t+dt [l] = l + nt[c*](DDl)- 1 dm t 
n[J]{Drf)- l dmt 
1 



cxp 



--n t [J]{Drf)-\t[c]dt 



(13) 
(14) 



Thus, 



^+dt[l] _1 = exv[-M<J}{DD^)- l dm t 

+ i 7 r t [c t ]( J D^ t )- 1 7r f [c]dt 

^l-T H [^]{DD^)- 1 dy t 

+ Tt t [c ] ]{DD ] )- 1 Tt t [c]dt (15) 

The infinitesimal evolution of a normalized conditional 
expectation is defined as 



ftt+dt[f] 

7Tt+dt[l] ' 



(16) 



Consequently, the evolution of the normalized condi- 
tional expectation is given by 

dn+dt[f] =*t[5f\dt + [7r t [(Vf)BD^ + / c t] 

- 7r t [/]7r t [c+]] (DD^-^dmt - w t [c}dt). (17) 

B. Linear filter 

Let us consider the conditional expectation of observ- 
ables at time t when the measurement outcomes M* are 
given. This is called filtering. The filtering equation is 
obtained by setting / = f(xt) in (|17| . Here we assume 
that the system is linear 



dxt = Axtdt + Bdivt, 
dmt — Cxtdt + Ddwt, 



(18a) 
(18b) 
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where A, B, C and D are constant matrices, and the ini- 
tial state Xo has a Gaussian distribution. 

By setting / — x t , we obtain the filtering equation 

d n [x t ] =ATT t [x t ]dt + (5C f + BD^(DD^)- 1 

x(dmt-Cir t [x t ]dt), (19) 

where 

S = Tr t [x t xl] - nt[x t ]nt[x]]. (20) 

For a Gaussian distribution, this is equivalent to 

5 = E[( ! r t -7rt[it])(a: t -7r t [!r t ])t], (21) 

so that S represents the error covariance matrix. Using 
(j!7p again, one can obtain the evolution of P as 

S =AS + SA* + BB^ 

- {SC ] + BD^)(DD^)~ 1 (SC^ + BD+)t. (22) 

C. Linear smoother 

Let us consider the conditional expectation of observ- 
ables at the initial time when the measurement outcomes 
M* are given. This is called smoothing (or fixed time 
smoothing in a strict sense.) The smoothing equation is 
obtained by setting / = f(xo) in (fT7|) . 

Here we assume the same linear system (|18|) . Setting 
/ = xq in (|17[) . we obtain the smoothing equation 

dTT t [x Q ] = KC^DD^idmt - Cn[x t ]dt), (23) 

where 

K =-K t [x x[] - 7r t {x ]irt [x\]. (24) 
For a Gaussian distribution, this is equivalent to 

K = E[(x - n t [x ]){xt - nt[xt]?]. (25) 
Using (|17p again, one can obtain the evolution of K as 

K =KA ] - /vC t (£'L' t )~ 1 [5C t + BD^]\ (26) 

where S is the covariance matrix of the filter (|20l . 

It is worth noting that smoothing has a different struc- 
ture from filtering because the filtering equation is self- 
consistent in the sense that it is updated by the current 
measurement outcome m t and conditional expectation 
7Ti[xi], whereas the smoothing equation depends on the 
filtering result. Due to this structural difference, while S 
is the error covariance matrix for filtering, K is not for 
the smoothing. In fact, K represents a correlation be- 
tween the current observable Xt and past one xq. It will 
be seen that K has a different meaning in a quantum 
case. 

For smoothing, the error covariance is defined as 

R = E[(x - ir t [x })(x - ^[a-o]) f ] (27) 



To calculate the time evolution of the error covariance, 
let us rewrite ([23]) as 

dn t [xo} = KC ] \DD^)- l [C{x t - n t [xt])dt + Ddw t ] (28) 

Then, it can be easily shown that 

R = -KC\DD J< )~ 1 CK Js . (29) 

III. QUANTUM FILTERING AND 
SMOOTHING 

The derivation of filtering and smoothing in the previ- 
ous section is different from a standard method which is 
based on the notion of projection and Wiener- Hopf equa- 
tion. We started from the conditional expectation for a 
nonlinear system. Then, filtering and smoothing equa- 
tions were easily obtained from the general framework 
in a unified way. In the quantum case, however, they 
should be considered separately because of noncommu- 
tativity. We start a general formulation from quantum 
filtering and consider a relation between the nondemoli- 
tion condition and smoothing. 

A. The system 

Let us consider two independent fields represented by 
symmetric Fock spaces Tjj (k = A, B) which are con- 
tinuously dilated spaces in t. We define field operators 
as 

C* = C k xm G L{T\), (k = A, B) (30) 

where L denotes a set of linear operators. These satisfy 
the quantum Ito rule 

dC k t dC l ? = 5 k idt, dC?dC l t = dC* ] ' dC\ = (31) 

which corresponds to In fact, a single quadrature of 
the field operator, e.g., the real part 

mf =Cf + (32) 

behaves as a classical Wiener process. The important 
property of the Wiener process, independency of incre- 
ments, follows from the structure of continuous dilation. 
The field operator at time t + dt is represented as 

Ct+dt =( ^x[o,t] ® C x [t-t+dt]- (33) 

Hence, future increments dCj? = C^ +dt — C t fe are indepen- 
dent of past field operators as 

lC^dC*] = [C^dC?]=0. for s<t. (34) 

In this sense, t is a dilation parameter rather than time 
for the fields T%. 
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Let us consider two quantum systems represented by 
Hilbert spaces Ji A and 3~C B which belong to Alice and 
Bob, respectively. Suppose that each system interacts 
with the field (k = A, B) independently. The total 
system is described as 



Jit 



(35) 



The dynamics of the system is described by a unitary 
operator Ut in this space, i.e., 

U t S L(K t ), (36) 

The time evolution of an arbitrary system operator 

Z AB e L(X A ® K B ) (37) 

is given by 



z AB = u\z AB 



U t eL(% t ). (38) 



Let 



L* e HCH k ) (39) 

be a system operator which couples to the field operators. 
Here we consider a unitary operator of the form 



dU t 



dC}L t 



L\dC t 



i 



L\L t dt 



where 



Lf 



Ct 



C A 
C B 



(40) 



(41) 



Here we ignore system Hamiltonians for simplicity. Since 
the two fields are independent, the unitary operator can 
be written as 



AB 



u A ® uf, 



(42) 

where U A and U B are independently defined in the same 
way as (|4T))) for each system. 

While the field operator C t fe can be thought of as a 
stochastic input to the system, the output of the system 
is defined by the field operator after the interaction with 
the system as 



D 



u}c k u t . 



(k = A, B) 



(43) 



Suppose that Alice makes a measurement of a single 
quadrature (real part) of her output field operator. Then, 
the measurement observable is represented as 



m t =D A + D 



At 



(44) 



Let W A and Mt be commutative von Neuman algebras 
generated by {u>^|0 < s < t} and {m t \0 < s < t}, re- 



spectively. From (|43|) , they are related to each other as 
M f = U}W A U t . (45) 
And also, any observable of Bob's system 

Z B etl{M B ) (46) 
is compatible with the measurement outcomes since 



[Z B ,M t 



for t> 0. 



(47) 



B. Nondemolition condition 

Unlike the classical case, the quantum conditional ex- 
pectation is not always well-defined. First of all, the 
measurement outcomes {m s |0 < s < i] can be defined as 
classically readable time-series data only when they are 
commutative with each other. This condition is naturally 
satisfied since for all s and t, 



[m s ,m t ] = U^ x{st) [w s ,wt}U ma ^ {S}t) = 0. 



(48) 



Secondly, for the conditional expectation of an observable 
/ to be well-defined as an estimate obtained from classical 
data M t , it should also be commutative with all M t , i.e., 



[/,M t 



0. 



This is called a nondemolition condition. From 
this condition is satisfied if / is of the form 



/ 
/ 



tAB 



t < Vr 
Vs < t < Vr 



(49) 

HZ}, 

(50a) 
(50b) 



The first case for r = t gives filtering and the second one 
corresponds to smoothing. We will consider these two 
cases separately in the following two subsections. (Here 
we do not consider prediction r > t because the result of 
prediction is given by a trivial master equation.) 



C. Quantum filtering 



Let us consider the first case / = Z AB of (|50[) . The 
quantum state equation corresponding to the classical 
one (fT]) is obtained by expanding the unitary operator in 
The resulting equation is written as 



df=3fdt+[L\dCt-dClL u fl 



dm t =(Lf + L^)dt + dw 



where 



If = L\fL t - l -L\L t f - \fL\L t 



(51a) 
(51b) 

(52) 



Our purpose is to calculate the conditional expectation 

7r t [/] = V[Z AB \M t ], (53) 

where P represents an expectation with respect to the 
initial density matrix of the system. By definition, we 
have 7Tf[/] g M t . 

To calculate the evolution of the conditional expecta- 
tion, we introduce a measure Q as 



P[Z AB ] 
where [16| 



Q[(V A ®U B yz AB (V A ®U B )}, 



dV/ 



L A dw A - L -L?L A dt 



Vt 



(54) 



(55) 
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Note that 



Q[(V t A ® U B )^Z AB (V A <8> f/f )|Wf ] e W 



(56) 



Since M t = U}WfU u we have 

p[^ B |M t ] = u}q[(v a ® u t B ^z^(y t A ® t/ t B )|w^, 

(57) 

which corresponds to ([?]) ■ Thus, the infinitesimal evolu- 
tion of the conditional evolution 7Tt [/] is obtained by cal- 
culating d[(V t A <g> U t B Yf(V t A ® C/f)] . Expanding V/ 1 ® C/ t B 
in the equation above, we have 

n t+dt [f] =n t [f] + ir t [?f]dt + n t [Lff + fL A ]dm t , (58) 

where TTt+dt denotes an unnormalized conditional expec- 
tation at time t + dt. 

Normalization requires the evolution of / = 1 as in the 
classical case and the resulting evolution of the normal- 
ized conditional expectation is given by 

ck t \f\ =nl?f]dt 



ir t [L A \f + fL A ]-MfML? + L?i] 



Ah 



x (dm t ~ ir t [L A + Lf]di) 



D. Quantum smoothing 



(59) 



cf 



Ci 



HD 







filter 









TTt[X t 



smoother 



-MX* 



FIG. 1: A schematic representation of quantum smoothing. 
Two systems, A and B, are initially entangled with each other, 
as expressed by ||. The output of Alice's system D A is de- 
tected by homodyne measurement, denoted by HD, which 
produces measurement outcomes m t . The filter produces the 
estimate 7r t [X t ]. The state of the filter is updated by the 
measurement outcome mt and its output 7Tt[Xt], as described 
by (1591 176[l . These two quantities also update the smoother, 
as shown in (|64l I82|) , and the output of the smoother is what 
we want, nt[Xo\- 



IV. QUANTUM LINEAR SYSTEM 

Before considering quantum smoothing in detail, we 
introduce a description of a quantum linear system in 
this section. Then, a smoothing problem for the linear 
system will be solved. 



Let us consider the conditional expectation for the sec- 

Unlike the classical case, 



ond case / = Z A <g> of flSj 
quantum smoothing cannot be obtained from (|59|) . How- 
ever, a modification for smoothing is rather simple. Our 
purpose here is to calculate 

P[Z A (g>Z B \M t 



n t[f] = r [ Zj t l JVL *J 
In this case, we introduce a measure as 



Q[v t 



M z A v A > 



(60) 
(61) 



Then, in a manner similar to the previous subsection, we 
can define the conditional expectation as 



n t [f] ~ UrQWrZfV* ® Z B \W A ]U t A (62) 

The infinitesimal evolution of the conditional expectation 
is given by 



dfr, [Z A 



-MHz? 



Z?]dt 



(63) 



+ n[(Lfz A + Z A L A )®ZZ]dm t 
Thus, the normalized conditional expectation obeys 
dir t [Zf®Z B ] =n t [?{Z A ). 



Z B ]dt 



M(Lfzt 



Z A L A ) ® Z B 



Mz? 



Z B ]n t [L A + Li 



Ah 



x (dm t - TTt [Lf + Lt J ]dt) . (64) 

Note that this is independent of Bob's dynamics. Once 
Alice starts measurement, her smoothing is not influ- 
enced by Bob's operation. An experimental scheme of 
smoothing is shown in FigfTJ 



A. Linear dynamics 

Assume that the two systems "K k (k — A,B) are 
bosonic and described by creation and annihilation oper- 
ators 

[af,a t fet ] = 1. (k = A,B) (65) 
Let us define orthonormal quadrature operators as 



Vi 



1 1 

— i i 



(66) 



If the operator L k in (|4H|) is linear in a k and et t , then 
we call the quantum system linear. Suppose that L\ is 
of the form 



Lt = 



if3 k 



7 fc + id k k 

o Ut i 



and the field operators are 



m\ 





i r 




■ c k - 




— i i 




. c? . 




1 1 " 




' D k ' 




— i i 







= m k . 



Then, a quantum linear system is expressed as 

x\dt 
dw k 



dx k 




' -A k /2 


-A k (G k )- 1 ' 


dm\ 




G k 


I 



(67) 

(68) 
(69) 

(70) 
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where 



G 



a k 7 fe 
(3 k S k 



A k = det G k 



(71) 



Let us assume that A fe 7^ and introduce non- 
orthogonal quadrature operators X k as 



X k = G k x\. 



(72) 



In this basis, the quantum linear system has a simple 
expression as 



-A k /2 



dX k 
dm k 



If we make a measurement of m\ = D k 
should be expressed as 



X k dt 
dw k 



fct 



dXl 
dm* 



' -A k /2 




~A k 


[[1 


oj 


I 


1 OjJ 


- A k 


B k 




' X k dt 


[c k 


D k 




dw k 



X k dt 
dw k 



(73) 
the system 

(74a) 
(74b) 



B. Quantum linear smoothing 

As in the classical case, quantum smoothing also re- 
quires the result of filtering, so we consider filtering for 
the quantum linear system first. Assume that the initial 
state of the system is an entangled Gaussian state. Let 
us define an operator-valued vector X t as 



X t = 



X? 



(75) 



From (|59|) . the expectation of X t conditioned on the mea- 
surement outcomes Mt is given by 



dn t [X t ] =An t [X t ]dt 
Here we have defined 

a AA 



F{dm t -CTr t [X t ]dt). (76) 



B 
D 



B J 



D A 
-1 



B B 
01, 



A B 

C=[C A ] . 
F = {SC^ - BD^){DD l[ ) 

and S is the error covariance matrix defined as 



S = 7T t [X t Xj] - TT t [X t ]7T t [Xj} 



s c s 1 



which obeys 

S =AS + SA^ + BB^ - F{DD^)F^ 



(77) 

(78) 
(79) 



(80) 



(81) 



In the quantum case, the smoothing equation is given 
by (|64p . For the linear system, the expectation of obscrv- 
ables X B conditioned on the measurement outcomes Mt 
is expressed as 

dir t [X B ] = K B C B ^ (DD^y 1 \dm t - C-K t [X t ]dt\ , (82) 



where the smoothing gain 

K B = n [X B X^] - n t [X B ]n t [X^] (83) 

satisfies 

K B = K B (A A ^ - C B \S A C A ^ + B A D A ^), (84a) 
K B (0) = S c (0). (84b) 

The error covariance matrix for quantum linear 
smoothing is defined as 



i?^E[efeft] 



(85) 



where ef = X B — TVt[X B ] is the smoothing error. The 
evolution of the error covariance is also obtained from 
dill) as 

R B = —K B C B \DD' { )~ 1 C B K B \ (86a) 
R B (0) = S B {0). (86b) 

V. EXAMPLE 



In this section, we consider an example of quantum lin- 
ear smoothing and optimization. For simplicity, assume 
that A = A A > and G B = I, i.e., the observable X B is 
equivalent to the orthonormal quadrature xf for Bob's 
system. (It can be seen that the case of A < has the 
same solution as A > 0. And A = will be discussed 
later.) 



A. Error covariance 



Let us express the error covariance matrix of filtering 



as 





S12 


Sl3 


S14 




S22 


023 


S24 


Sis 


S23 


S33 


S34 


. Su 


S24 


S34 


S44 . 



s = 



From (|81j). Su obeys 

Su = -AS 11 + A 2 -(S 11 -A) 2 . 
Its solution is given by 

A 



(87) 



Su = 



1 — /ue 



-At ' 



where 



From 



H = 1 - 



Sii(0) 

the smoothing gain obeys 



K B = - ^K B -K B 



Su — A S12 




(89) 



(90) 



(91) 
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The error covariance of filtering (|89|) leads to solutions 



= 5i 3 (0)e-T' 
Jf£ = S u (0)e-T' 



1-M 
i - ^e" A * : 

i -^e- At ' 



(92a) 
(92b) 



We are particularly interested in the smoothing error 
covariance matrix here because it is an important quan- 
tity to evaluate the performance of smoothing. From 
we have 



R B (t) =S B {0) 



where 



S? 3 (0) Si 3 (0)5i 4 (0) 
5i 4 (0)5 13 (0) S? 4 (0) 



h(t) 



1 - e 



-A/ 



1 — /ie 



-At ' 



s u (oy 

(93) 



(94) 



B. Measurement optimization 

The reduction in the smoothing error is given by the 
second term of (|93|) . This is obviously related to the 
amount of information about Bob's initial state that 
Alice can obtain through the continuous measurement. 
Thus, it is natural to maximize this quantity by choosing 
measurement parameters. To this end, let us introduce 
I(t) as 



/(*) 



-h(t) 



Sii(O) 

where for a non-negative constant matrix 



e = [ 5i 3 (o) s u (o) ] e 



5i 3 (0) 
5i 4 (0) 



(95) 



(96) 



Since I(t) is monotonic in time, the efficiency of esti- 
mation in an early stage is determined by 



<U_ 



(97) 



On the other hand, I(t) approaches asymptotically to 



It should be noted that G A is determined by two factors: 
What Alice measures and the coupling constant between 
Alice's system and the external field. Here we assume 
that the initial entangled state is a two mode squeezed 
state and investigate the smoothing performance in de- 
tail. For the purpose of smoothing, it is natural that she 
wants to minimize Fi[(x B —n t [x B ]) 2 ] and E[(y B — Kt[Vo ]) 2 ] 
at the same time. This corresponds to taking a diagonal 
weight matrix 9 = diag[9i 62] in ([96]) . 

In the orthonormal basis (note that G B = I now) 



X 



(99) 



the covariance matrix of the two mode squeezed state is 
given by 



cosh(r)/ 2 x2 sinh(r)J 
sinh(r)J cosh(r)J 2X 2 

where r is a squeezing parameter and 

-1 



J 



1 



(100) 



(101) 



Note that the covariance matrix is invariant under local 
shift operations. Now, since only Alice's system is in the 
non-orthogonal basis X A , the initial covariance matrix 
S(0) is given by 



S(0) 



cosh(r)G A G At sinh(r)G A J 
sinh(r)JG' 4t cosh(r)/ 2 x2 



(102) 



Thus, the two quantities representing the smoothing per- 
formance are calculated as 



dl 



[6 1 (a A ) 2 + e 2 ( 7 A ) 2 ]sinh 2 (r), 



t=o 



J(oo) = 



e 1 (a A ) 2 + 9 2 ( 7 - 4 ) 2 sinh"(r) 
(a A ) 2 + (7 A ) 2 cosh(r) 



(103) 
(104) 



Note that a A and ^y A are related to Alice's measurement 
observable because L A + L Aj( = a A x A +j A y A . (This also 
includes a coupling constant between Alice's system and 
the field implicitly.) However, if she wants to equally 
estimate x B and y B , then the asymptotic performance 
depends only on the initial squeezing. 



7(oo) 



Sn(0)' 



(98) 



A good smoother or good communication between the 
two parties is therefore defined by these two quantities. 
Alice will design her measurement to maximize the first 
one if she wants to estimate Bob's initial state as soon as 
possible, whereas the second one should be maximized if 
she wants a precise estimate after a long period of mea- 
surement. 

These two quantities depend on the matrix G A and 
the initial entanglement shared between the two parties. 



C. Smoothing uncertainty relation 

Since the system is rather simple in the present case, 
the smoothing performance can be directly seen by the 
smoothing error covariance matrix R B (t) . From (|93[) , we 
have 



Rf^t) = cosh(r) 



Rg(t) =cosh(r) 



{a A Yh(t) sinh>) 
(a A ) 2 + (7 A ) 2 cosh(r) 

j-f A ) 2 h(t) sinh 2 (r) 
(a A ) 2 + (7 A ) 2 cosh(r) 



(105a) 
(105b) 
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Thus, we can reduce the smoothing error of x B and y B 
to zero simultaneously by preparing perfect two mode 
squeezing r — > oo and a A ^ 0, j A ^ 0. And also, from 
(|93p again, we have the smoothing uncertainty relation 



det R B (t) = cosh 2 (r) - h(t) sinh 2 (r) > 1, 



(106) 



where the equality is attained when t — > oo. 

The mutual information of x B and m f is given by 



I[x*;m t ] = -]og\S B (0)(R*(t))- 1 \ 



1 



log(l - h(t) tanh 2 (r)). 



(107) 
(108) 



It can be decomposed into 

I[x B ;m t ] =llog|5 B (0)(S' B (t))- 1 | 
+ l\og\S B (t)(R B (t))- 



(109) 



The first term represents the amount of information 
about x B obtained from the filter, whereas the second 
term is an information difference between the filter and 
smoother. This decomposition actually clarifies a differ- 
ence between the roles of filtering and smoothing in FigJT] 
For a weak squeezing limit r« 1, we have 



I[x B ;m t ] 



'— log — 

2 8 1 



1 



1, 1- 
- log 

2 B 1 



"*/i(t)r 2 



h{t)i 



(110) 



Note that we have set G B = I. The first term decays by 
the interaction between Bob's system and the external 
field at the rate of e - ( dotG Then, the information 
provided by the smoother complemcntarily increases, as 
shown by the second term. By contrast, in the early stage 
of smoothing t <C r with a strong squeezing r » 1, the 
mutual information is expressed as 



I[x B ;m t ) 



t-iio g (i-M*)) 



- 1. 



(in) 



If we have entanglement strong enough, a simple filtering 
P[xf |M t ] gives a good estimate of x B , as shown in the 
first term, and the extra structure (the smoother) is not 
necessary. However, after a certain period of time, the 
smoother starts working normally and the filter does not 
produce information on x B . 



VI. DISCUSSION 

So far, we have considered the case of A > 0. It is not 
difficult to see that the same results hold true for A < 0. 
However, A = is a singular situation and we need to 
investigate this case separately. To see the singularity, 
we assume that G A is of the form 



G' 



a A a A 

(3 A (3 A 



(112) 



Then, 



Xi 



G' 



x? 
Vt 



a A (x A + y A ) 
A (x A +y A ) 



(113) 



which indicates that the two non-orthogonal quadratures 
are degenerate. Hence, the system is expressed by the 
orthonormal basis x A as in (l70"1). 



(114) 



Let us introduce a similarity transform 



x 



Al 



Txi 



T = —= 



1 

3 



1 1 

1 -1 



Since T is unitary, the commutation relation is invariant 



under the transform, i.e., [x 
system is now expressed as 



Vt 



dxf = 4= 

V2 
1 

3 



dm? = — 





-2b A 

2a A 
2b A 





2a A 



x?'dt 



dwi 



Alice's 

(115a) 
(115b) 



The first equation implies that x?' (the first element of 
x?') is static, i.e., x?' is back-action free. Yet, Alice can 
obtain information about x?' by measuring m? (the first 
element of mf.) This is quantum nondemolition (QND) 



measurement of x?' 



x A 



y? . In fact, it follows from 



the unitary operator (|40|) that 



.At 



is compatible with 



Thus, in this case, Alice can perform smoothing on 



her observable x?' . In other words, nt[x A/ ] = P[x A |Mf] 
is well-defined. This is, however, almost trivial because 
xf' = x Al . The filtering equation for ^[a;^'] gives the 
smoothing result ir t [x Al ] . 



VII. CONCLUSION 

We have introduced a general formulation of quantum 
smoothing. The idea of smoothing cannot be defined 
for quantum systems in a straightforward manner as in 
the classical case because of the nondemolition condition. 
Quantum smoothing is possible only through entangle- 
ment except for QND measurement discussed above. The 
initially shared entanglement is used for the estimation 
of the initial state, and the state after smoothing asymp- 
totically becomes completely separable. In this sense, 
the mutual information is related to the amount of en- 
tanglement in the initial state. It is worth noting that 
the mutual information given here is independent of the 
detailed structure of the matrix G k . On the other hand, 
other measures of entanglement such as the negativity are 
sensitive to Alice's estimation method as the smoothing 
parformance was dependent on a A and 7 . 



9 



[1] L. Thomsen, S. Mancini, and H. M. Wiseman, J. Phys. 

B 35, 4937 (2002). 
[2] J. Geremia, J. K. Stockton, and H. Mabuchi, Science 

304, 270 (2004). 
[3] J. Geremia, Phys. Rev. Lett. 97, 073601 (2006). 
[4] A. Negretti, U. V. Poulsen, and K. Molmer (2007), quant- 

ph/0703040. 

[5] N. Wiener, Extrapolation, Interpolation, and Smoothing 
of Stationary Time Series (Wiley, New York, 1949) . 

[6] R. E. Kalman, Transactions of the ASME-Journal of Ba- 
sic Engineering 82, 35 (1960). 

[7] C. W. Gardiner and P. Zoller, Quantum Noise (Springer, 
Heidelberg, 2004). 

[8] M. Yanagisawa and H. Kimura, IEEE Trans. Autom. 
Contr. 48, 2107 (2003). 

[9] E. B. Davies, Quantum Thoery of Open Systems (Aca- 



demic Press, London New York San Francisco, 1976). 
[10] V. P. Belavkin, J. Multivar. Anal. 42, 171 (1992). 
[11] L. Bouten, R. van Handel, and M. James (2006), quant- 

ph/0703040. 

[12] H. J. Carmichael, An Open Systems Approach to Quan- 
tum Optics (Springer- Verlag, Berlin Heidelberg New 
York, 1993). 

[13] A. Griffin, D. W. Snoke, and S. Stringari, eds., Bose Ein- 
stein Condensation (Cambridge University Press, Cam- 
bridge, 1995). 

[14] M. W. Jack, M. J. Collett, and D. F. Walls, Phys. Rev. 

A 54, R4625 (1996). 
[15] J. E. Lye, J. J. Hope, and J. D. Close, Phys. Rev. A 67, 

043609 (2003). 
[16] A. Holevo, J. Soviet Math 56, 2609 (1991). 



